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ABSTRACT 

The  concept,  instrumentation,  and  analytical  techniques  for  the 
investigation  of  attenuation  of  wind  generated  water  surface  waves  in 
a  small  scale  real  environment  by  separation  of  generation  and  propaga- 
tion by  time  in  the  dying  wind  situation  were  developed,  and  an  empirical 
attenuation  function  applicable  to  a  wide  range  of  frequencies  proposed. 
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I.   INTRODUCTION 

The  concept  of  wind-driven  ocean  waves  as  spectra  of  harmonic 
components  which  follow  the  laws  of  linear  hydrodynamics  was  introduced 
in  1948  by  Barber  and  Ursell  (1).   There  has  since  been  considerable 
progress  toward  the  mathematical  and  physical  description  of  ocean  wave 
processes,  with  the  ultimate  objective  being  the  prediction  with  reason- 
able certainty  of  future  wave  conditions  at  any  geographical  point  on 
the  basis  of  worldwide  synoptic  observations  of  both  wave  and  wind  fields 
A  summary  of  the  major  contributions  to  knowledge  of  ocean  waves  was 
compiled  in  1967  by  Cartwright  (2). 

While  theoretical  models  of  wave  generation,  dispersion  and  geo- 
metrical spreading  have  been  fairly  well  confirmed  empirically  and 
are  used  as  the  current  basis  of  wave  forecasting,  the  attenuation  of 
waves  during  propagation  is  less  well  understood;  and  an  attenuation 
function  is  consequently  not  applied  in  many  wave  —forecasting  techni- 
ques. 

Some  classical  theories  of  wave  attenuation,  namely  air  resistance 
and  turbulent  friction,  were  critically  examined  by  Darbyshire  (4)  and 
found  wanting  as  singular  explanations.   The  theory  of  non- linear  scat- 
tering proposed  by  Phillips  (6)  appears  promising  and  has  been  reason- 
ably supported  by  Snodgrass,  et.  al.  (11) .   Conceivably,  the  answer 
may  lie  in  a  combination  of  these  processes. 

As  has  been  the  case  in  studies  of  other  environmental  processes, 
it  appears  that  the  phenomenon  under  consideration  must  be  completely 
described  empirically  before  the  question  of  theory  can  be  resolved. 
So  far,  environmental  studies  of  attenuation  have  been  confined  to  low 
frequency  ocean  swell,  where  the  range  of  measurable  frequencies  is 
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included  in  only  one  of  several  logarithmic  decades  of  possible  surface 
wave  frequencies.   The  classical  theories  and  the  generalized  model 
developed  later  in  this  paper  call  for  the  attenuation  function  to  be 
exponential  in  propagation  distance  and  frequency  dependent.   When  the 
generalized  model  is  applied  to  ocean  wave  statistics,  the  data  are 
scattered  within  such  a  small  band  of  the  frequency  scale  that  func- 
tions supporting  almost  all  proposed  theories  can  be  fitted  equally  well 
to  the  data. 

Accordingly,  the  primary  objective  of  this  research  was  the  collec- 
tion of  attenuation  data  in  a  small  scale  real  environment  where  the 
measurable  frequencies  are  considerably  displaced  from  those  normally 
measured  in  ocean  swell  in  order  to  give  better  definition  to  the  at- 
tenuation function.   This  requirement  for  collection  of  data  gave  rise 
to  the  second  objective,  the  development  and  evaluation  of  a  system  for 
recording  and  analyzing  waves  in  a  small  scale  real  environment. 
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II.   CONCEPTS 

A.    COMPARISON  OF  ENVIRONMENTS 

Ocean  wave  attenuation  studies  involve  the  measurement  of  the  energy 
of  the  relatively  low  frequency  components  of  swell  propagated  over  long 
distances  from  the  area  of  generation.   It  is  known  that  these  low 
frequencies  required  high  winds  of  appreciable  duration  for  their  develop- 
ment, hence  the  source  of  the  swell  waves  can  be  located  with  fair  ac- 
curacy in  relation  to  the  distance  from  the  measuring  instrument.   A 
knowledge  of  wind  or  pressure  conditions  in  the  fetch  can  provide  a 
reasonable  estimate  of  the  energy  of  waves  in  the  fetch. 

High  frequency  waves  can  be  generated  by  the  prevailing  low  veloc- 
ity winds  which  are  normally  found  over  most  of  the  ocean  surface.   For 
this  and  other  reasons  of  distance  and  instrumentation,  it  would  be  ex- 
tremely difficult  to  make  meaningful  measurements  of  these  frequencies 
with  sufficient  precision  in  the  ocean  to  determine  attenuation.   For 
purposes  of  this  paper,  high  frequencies  shall  be  defined  as  those  fre- 
quencies approaching,  but  not  including,  frequencies  in  which  surface 
tension  is  the  predominant  restoring  force  rather  than  gravity. 

A  lake  or  non-tidal  estuary  of  moderate  depth  provides  an  ideal 
facility  where  experimental  control  is  possible  in  the  comparative 
measurement  of  high  frequency  deep  water  waves,  although  the  scale  of 
wind  circulation  in  relation  to  the  size  of  the  lake  or  estuary  provides 
no  geographical  separation  between  areas  of  generation  and  non-genera- 
tive propagation.   On  this  scale,  however,  the  mean  wind  may  change 
rapidly  in  speed  and/or  direction.   In  a  case  where  the  wind  speed 
diminishes  rapidly,  some  frequencies  generated  at  the  higher  wind  speed 
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would  receive  no  more  energy  and  would  propagate  as  swell.   This  con- 
cept, that  fetch  and  swell  can  be  separated  by  time  just  as  there  is 
the  geographical  separation  in  the  ocean,  is  a  fundamental  assumption 
of  this  paper. 

Bo    THE  GENERALIZED  ATTENUATION  FUNCTION 

Let  S(f,x)  represent  the  energy  density  spectrum  computed  from  a 
wave  record  at  a  downwind  (propagation)  position,  and  let  S(f,0)  represent 
the  energy  density  spectrum  computed  from  a  wave  record  at  an  upwind 
(generation)  position,  where  f  is  the  component  frequency  and  x  is  the 
separation  between  two  points  of  measurement  along  the  axis  of  maximum 
propagation  (the  direction  of  the  mean  wind  during  generation).   Then, 

S(f,x)  -  S(f,0)   k(x,0,f)   G(f,x)  (1) 

where  G(f,x)  is  the  attenuation  function,  and  k(x,9,f)  denotes  what 
portion  of  S(f,0)  would  reach  the  downwind  position  if  no  attenuation 
were  involved.   This  is  a  function  of  geometry  and  the  beam  pattern  in 
the  generating  area.   In  this  paper  the  model  of  Pierson,  Neumann  and 
James  (9)  will  be  used,  where  for  all  f 

Q, 


/  4    2 
k(x,0,f)  =  — ~-       COS  9  d9  (2) 

77     S 


The  angles  9  and  9.  are  measured  at  the  edges  of  the  fetch  from  the 
direction  of  the  axis  described  above  to  the  direction  of  the  downwind 
measurement  point. 

Let  P(f,x)  be  a  prediction  of  the  spectrum  at  the  downwind  position 
which  has  not  allowed  for  attenuation.   Then, 

P(f,x)  =  S(f,0)   k(x,9,f)   H(f,x)  (3) 
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where  H(f,x)  is  an  error  function  accounting  for  possible  statistical 
errors  in  measurement  at  the  upwind  and  downwind  positions,  and  in- 
cludes differences  between  the  angular  spreading  model  and  the  actual 
angular  spreading  phenomenon.   The  model  has  proved  to  be  fairly  representa- 
tive of  actual  geometrical  spreading,  and  H(f,x)  should  be  at  most  a  very 
slowly  changing  function  of  x.   The  spacial  variation  of  H  is  thus  ignored. 
Divide  (3)  by  (1) 

Rff  x\    =  -ULiSl  =   H(f,x) 

R(f'X)    S(f,x)     G(f,x)  (4) 

Differentiate 

R'(f,x)  «— £_  =   -H(f,x)   G'(f'x)  (5) 

d  G^(f,x) 

Substitute    (4)    in   (5) 

R'(f.x)    =      -R(f,x)        [j^ffi  (6) 

Now,  assume  that  R(f,x)  can  be  represented  by: 

R(f,x)  =  C  emX  (7) 

where  m  =  m(f),  and  C  is  a  constant.   This  assumption  appears  to  be 
valid  from  empirical  data  in  ocean  swell  studies  over  a  wide  range  of  x. 
Differentiating  (7)  and  substituting  (7)  in  the  result, 

R'(f,x)   =  Cm  emX  =  m  R(f,x)  (8) 

then  substitute  (8)  in  (6),  and 

G'(f,x)  =   -m  G(f,x)  (9) 


17 


which  is  a  simple  differential  equation  with  solution 

G(f,x)   =  GQ(f,0)  e"mX  (10) 

Clearly,  at  the  upwind  measuring  position,  the  function  G  (f,0)  =  1 


and 


G(f,x)  =  e~mX  (11) 


The  attenuation  coefficient  m  is  a  function  of  frequency  and  environ- 
mental conditions.   Empirical  evidence  indicates  that  it  increases  with 
frequency,  and  from  physical  considerations  it  should  reasonably  in- 
crease from  following  wind  to  opposing  wind  conditions.   Other  environ- 
mental conditions  may  also  be  important. 
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III.   INSTRUMENTATION 

Two  resistance  wave  meters ,  adapted  from  a  design  by  Pearlman  (5), 
were  constructed  and  used  to  collect  data  for  this  paper.   A  schematic 
diagram  of  the  circuit  used  in  both  meters  is  shown  in  Figure  1. 

Pearlman 's  paper  demonstrated  that  this  type  of  wave  sensor  is  es- 
sentially linear  over  a  range  of  wave  heights  of  a  few  inches  in  static 
and  vertical  oscillation  tests  and  has  errors  of  less  than  one  percent 
when  tested  in  orbital  motion.   Static  tests  of  the  modified  design  used 
in  this  research  confirmed  the  linear  characteristics  of  the  instruments, 

The  fixed  resistors  in  the  wheatstone  bridge  were  standard  57o  toler- 
ance resistors ,  but  were  matched  by  measurement  on  a  resistance  bridge 
to  have  values  within  0.5  ohm  of  the  resistance  shown.   The  variable 
resistor  was  used  to  null  the  bridge  at  an  initial  immersion  in  the 
field  prior  to  recording  data. 

The  signal  generators  used  were  Hewlett  Packard  Model  2000  audio 
oscillators.   The  recorder  was  a  14  channel  Ampex  CP-100  tape  recorder. 
As  the  signal  generators  and  the  tape  recorder  both  have  one  grounded 
terminal,  the  isolation  transformer  was  required  to  prevent  grounding 
one  leg  of  the  wheatstone  bridge. 

Alternating  current  was  used  to  minimize  electrolytic  effects  on 
the  electrodes.   The  full  wave  rectifier  shown  in  Figure  1  demodulated 
this  carrier  frequency s  passing  a  direct  current  signal  to  the  re- 
corder that  varied  in  response  to  the  waves.   The  signal  generators 
were  operated  at  3  volts  (RMS)  at  about  2000  Hz;  however,  when  operat- 
ing two  wave  meters  on  one  recorder,  it  was  necessary  to  have  about  150 
Hz  difference  in  the  output  of  the  signal  generators  to  prevent  a  beat 
frequency  from  interfering  with  the  wave  records. 
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A  dual  beam  oscilloscope  was  used  to  monitor  the  input  and  output 
of  the  tape  recorder  to  insure  proper  recording  of  the  wave  signal.   The  1 
Mfd  capacitor  used  across  the  rectifier  output  matched  the  input  impedance 
of  the  tape  recorder  and  oscilloscope,  but  the  use  of  other  equipments 
with  different  input  impedances  would  require  a  change  in  capacitance  to 
provide  a  smooth  rectifier  output  without  attenuating  wave  frequencies. 

The  Ampex  CP-100  tape  recorder  was  operated  with  frequency  modula- 
tion recording  and  reproducing  amplifiers  for  the  wave  record  channels 
at  a  tape  drive  speed  of  1-7/8  inches  per  second.   A  third  channel  was 
used  with  amplitude  modulation  (direct)  amplifiers  for  a  voice  time 
record. 

Details  of  the  sensing  probes  are  shown  in  Figure  2.   The  frames 
were  designed  to  provide  adequate  rigidity,  yet  to  minimize  turbulence 
and  reflections.   The  lower  horizontal  supports  were  designed  to  be  at 
a  depth  greater  than  one-half  wavelength  of  the  anticipated  waves.   The 
fishing  line  shown  was  used  as  an  insulator  and  to  minimize  the  amount 
of  electrode  wire  immersed.   If  the  electrodes  had  extended  to  the  lower 
horizontal  support,  the  wheatstone  bridge  could  not  have  been  brought 
to  near-balance  by  the  variable  resistor.   In  practice,  it  was  found 
that  the  bridge  could  be  nulled  at  an  electrode  immersion  of  no  more 
than  2  centimeters. 

The  probe  frames  were  attached  to  1  inch  diameter  galvanized  pipes 
driven  into  the  lake  bottom  by  means  of  the  1\   inch  diameter  sleeves 
shown  in  Figure  2.   The  meter  circuits  were  mounted  on  top  of  the  frame 
and  connected  to  the  oscillators  and  recorder  ashore  with  Alpha  #2311 
low  capacitance  shielded  microphone  cable  laid  on  the  bottom. 
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IV.   DATA  COLLECTION 

Roberts  Lake,  1500  feet  from  Monterey  Bay  in  Seaside,  California, 
was  used  for  the  collection  of  data  in  support  of  this  paper.   The  lake, 
shown  in  Figure  3,  is  uniformly  5  to  5h   feet  deep,  with  a  soft  mud  bottom. 
Dense  reed  growth  in  the  water  around  the  circumference  minimizes  wave 
reflections. 

The  Monterey  Bay  area  was  considered  suitable  for  this  small  scale 
attenuation  study  because  its  climatological  characteristics  are  favor- 
able to  sudden  changes  in  the  wind.   In  addition  to  unpredictable  changes 
in  the  wind  field,  diurnal  heating  and  cooling  of  the  Salinas  Valley  20 
miles  inland  results  in  a  sea  breeze  during  the  day  which  dies  suddenly 
at  dusk.   This  general  characteristic  is  modified  or  not  apparent  when 
the  central  California  coast  is  under  the  influence  of  transient  high  or 
low  pressure  regions.   In  three  weeks  of  daily  attempts  to  collect  data, 
only  one  day  produced  the  desired  rapidly  dying  wind. 

Two  sets  of  data,  identified  as  data  sets  A  and  B,  were  collected 
on  21  November  1968.  The  wind  speeds,  as  measured  with  a  Casella  pre- 
cision anemometer  mounted  eight  feet  above  the  lake  surface,  are  shown 
in  Figures  4  and  5.  The  Casella  anemometer  provides  a  mean  wind  speed 
by  recording  the  revolutions  of  a  three  cup  rotor  on  a  digital  counter. 
The  wind  data  in  data  set  A,  Figure  49  were  recorded  for  one-minute  aver- 
ages, and  in  data  set  B,  Figure  5,  for  half-minute  averages. 

Details  of  wave  sensor  placement  are  shown  in  Figure  3.   For  data 
set  A,  the  sensing  probes  were  on  an  axis  of  230° (M) .   Although  the 
actual  distance  between  probes  was  400  feet,  the  direction  of  the  wind 
(from  300° (M))  reduced  the  effective  propagation  distance  to  138  feet. 
For  data  set  B,  the  probes  were  120  feet  apart  on  an  axis  of  290° (M). 
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Prior  to  collecting  the  wave  data,  both  probes  were  immersed 
simultaneously  in  relatively  still  water  near  the  edge  of  the  lake  at 
two  levels  of  immersion  and  one  minute  recordings  of  the  signals  were 
placed  on  magnetic  tape  for  calibration  of  the  wave  meters.   The  average 
values  of  the  one  minute  calibration  signals  were  later  used  to  scale 
the  wave  records  to  actual  amplitudes  in  millimeters. 
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V.      ANALYSIS   OF   DATA 

A.  DIGITIZING 

The  calibration  records  and  selected  portions  of  the  wave  records 
were  digitized  with  an  SDS  9300-CQMCQR  5000  hybrid  computer  system. 

The  Ampex  CP-100  tape  recorder  was  patched  to  the  COMCOR  5000  analog 
computer  through  operational  amplifiers  included  in  the  logic  board  of 
the  computer  and  a  low  pass  filter  constructed  from  logic  board  components. 
The  amplifiers  were  adjusted  to  amplify  the  maximum  signal  on  the  magnetic 
tape  record  to  slightly  less  than  the  maximum  input  voltage  of  the  analog 
computer. 

A  schematic  diagram  and  plot  of  the  measured  characteristics  of  the 
low  pass  filter  are  shown  in  Figure  6.   To  minimize  aliasing  in  later 
spectral  calculations,  128  Hz  was  selected  as  a  filter  cut-off  frequency. 
Although  frequencies  of  that  order  were  far  beyond  the  range  of  expected 
surface  wave  frequencies,  noise  could  possibly  have  existed  in  that  range 
from  either  inherent  equipment  noise  or  beats  between  the  wave  meter 
oscillators. 

To  provide  a  minimum  of  two  samples  per  cycle  up  to  the  cut-off 
frequency,  the  record  was  sampled  at  a  rate  of  256  samples  per  second. 
The  digital  computer  sampling  program  used  was  developed  by  Post  (10). 

B.  CONVERSION  OF  DATA  FOR  IBM  360  COMPUTER  SYSTEM  COMPATIBILITY 

The  core  capacity  of  the  SDS  9300  computer  used  for  digitizing  the 
data  was  far  too  small  for  spectral  analysis  of  the  wave  records  (see 
next  section).   Therefore,  use  was  made  of  the  much  larger  core  capacity 
of  the  available  IBM  360  computer  system.   The  two  computers  do  not  share 
common  symbology  and  storage  parameters;  however,  a  program  developed  by 
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Post  (10)  for  conversion  of  data  from  one  system  to  the  other  was  used 
to  transfer  and  convert  the  digitized  data  from  the  tape  record  made  on 
the  SDS  9300  computer  to  a  tape  record  compatible  with  the  IBM  360 
computer  system. 

C.    CALCULATION  OF  ENERGY  DENSITY  SPECTRA 

In  a  particular  finite  segment  of  a  wave  record  of  length  T  ,  let 
A(t,  )  represent  the  time  series  of  N  discrete  wave  amplitudes  which  are 

K 

spaced  at  intervals  of  time  At.   Then  the  time  of  each  sample  amplitude 
is 

tk  =  k  At,  k  -  1,...,N         (12) 

The  record  is  transformed  from  the  time  domain  to  the  frequency 

domain  by  the  Fourier  transform 

N-l 
1<V--1-    I    A(tk)   e-2TTiJk/N  ,j.  0.....N-1    (!3) 

k  =  0 


where   f .  =  j  Af  a*"e  the  component  frequencies 
A(f.)  are  the  complex  amplitudes  of  the  f . 

and      i  =      -fT~ ' 

In  the  computer,  this  transform  is  accomplished  by  the  IBM  library 
subroutine  HARM,  which  uses  the  fast  Fourier  transform  algorithm  developed 
by  Cooley  and  Tukey  (3) 9  and  requires  that  N  be  an  integral  power  of  2. 

A(f .)  is  a  two-sided  amplitude  spectrum  where  N/2  corresponds  to  the 

aliasing  or  Nyquist  folding  frequency,  f  ,  which  is  the  highest  frequency 

c 

analyzed.   If  the  wave  record  had  component  frequencies  above  f  ,  the  ampli' 
tudes  of  these  frequencies  would  be  aliased  into  the  spectrum  A(f.).   A 
Nyquist  frequency  of  128  Hz  was  selected  on  the  basis  of  the  analog  filter 
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characteristics  as  discussed  in  section  A.   The  Nyquist  frequency  deter- 
mines the  minimum  sampling  interval,  At,  by  the  relation 

At  =  l/2f  (14) 

c 

The  energy  density  spectrum,  S(f .)  is  given  by 

S(f.)   -   2  \kf.)Uts)\    >  ,.„ m  (15) 

A£ 

where  A*(f .)  is  the  complex  conjugate  of  A(f  .)  and  the  factor  2  accounts 
for  the  inclusion  of  only  half  of  the  two-sided  symmetrical  amplitude 
spectrum.   The  S(f.)  represent  the  average  energy  per  unit  frequency  in 
an  elementary  frequency  band  of  width  Af  centered  on  f . . 

The  S(f .)  of  one  realization  of  an  ergodic  random  process  is  a  set 
of  statistics,  each  of  which  is  considered  to  have  the  Chi-square  distri- 
bution with  two  degrees  of  freedom.   The  number  of  degrees  of  freedom 
determines  the  confidence  limits  for  a  statistic,  and  a  meaningful  compari 
son  of  statistics  requires  that  the  range  of  confidence  limits  be  suffi- 
ciently small.   The  relation  between  the  number  of  degrees  of  freedom 
and  the  confidence  limits  for  the  Chi-square  distribution  is  given  by 
Pierson  and  Marks  (8). 

An  accepted  method  of  increasing  the  confidence  in  the  statistical 
estimate  of  energy  density  is  to  average  a  number,  M,  of  elementary 
frequency  bands  into  a  wider  equivalent  band  of  width  B  .  When  this  is 
done,  each  bandwidth  B  has  2M  degrees  of  freedom,  As  the  two  degrees 
of  freedom  in  each  elementary  bandwidth  reflect  two  independent  samples 
in  the  time  domain,  the  record  length,  T  ,  required  for  the  analysis  is 
then  specified  by: 

T  =   N  At   =  2  f  M  At   /  B  (16) 
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Obviously,  the  increased  confidence  is  bought  at  the  expense  of  defini- 
tion,, 

The  parameters  used  for  analysis  of  the  wave  records  collected  on 
21  November  1968  are  summarized  as  follows; 

T  =  64  seconds 
r 

N  =  16384 

£t  =  1/256  seconds 

f   =  128  Hz 
c 

B  -  0.5  Hz 
e 

M  =  32 

For  64  degrees  of  freedom,  the  957»  confidence  limits  are  1.5  and 

0.7,  which  means  that  for  an  infinite  number  of  realizations  of  this 

particular  wave  process,  95%  of  the  realizations  of  the  spectrum  would 

fall  between  1.5  S(f  .)  and  0.7  S(f.). 

J  J 

In  the  discussion  of  the  results  of  the  spectral  analyses  in  the 
following  section,  confidence  limits  have  been  omitted  for  clarity,  but 
their  existence  must  be  recognized.   In  fact,  the  computed  S(f .)  are  the 
most  likely  values. 
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VI.   RESULTS  AND  IMPLICATIONS  OF  SPECTRAL  CALCULATIONS 

In  the  text  and  figures  of  the  discussion  that  follows,  the  time 
(hours :minutes: seconds)  assigned  to  a  spectrum  is  the  mid- time  of  the 
64-second  wave  record  from  which  the  spectrum  was  derived.   For  the  wind 
data  subset  of  data  sets  A  and  B,  refer  to  Figures  4  and  5,  respectively, 

A.  COMPARISON  OF  SIMULTANEOUS  UPWIND  AND  DOWNWIND  SPECTRA 

Figure  7  shows  the  energy  density  spectra  computed  from  simultan- 
eous realizations  at  the  upwind  and  downwind  wave  probe  positions  at 
11:55:56  and  11:57:00.   This  comparison  implies  that  spectra  from  down- 
wind probe  data  may  be  underestimated  by  a  factor  of  approximately  two. 
The  most  obvious  cause  for  this  under-estimation  is  an  error  in  calibra- 
tion* however,  it  is  possible  that  the  difference  may  be  an  effect  of 
the  geometry  of  the  lake. 

B.  ATTENUATION  INDICATED  BY  DATA  SET  A 

Figure  8  shows  the  energy  density  spectra  at  the  upwind  probe  posi- 
tion at  11:57:00  when  the  average  wind  at  the  eight  foot  level  had  been 
four  knots  or  more  for  at  least  30  minutes,  and  at  the  downwind  probe 
position  112,  176,  and  240  seconds  later.   During  the  four-minute  inter- 
val represented,  the  wind  diminished  to  1.3  knots. 

In  the  development  of  the  generalized  attenuation  function  earlier 
in  this  paper,  the  true  energy  density  spectra  at  the  generation  and 
propagation  positions  were  related  by; 

S(f,x)  =  S(f,0)   k(x,9,f)   G(f,x)  (1) 

and  the  attenuation  function  was  specified  by: 

G(f,x)  =  e"mX  (11) 


33 


UPW  NO 

DOV  NWIND 


55-56 


FREQUENCY,    SEC 


-I 


FIGURE    7.     COMPARISON     OF     SIMULTANEOUS        SPECTRA 


34 


.003- 


FREQUENCY.     SEC 


FIGURE    8.     ENERGY     DENSITY    SPECTRA- DATA    SET     A 

35 


If  it  is  assumed  that  the  energy  density  spectra,  S(f  .,0)  and 

J 

S(f  .,x),  derived  from  measurements  of  the  upwind  and  downwind  probe 
positions,  respectively,  are  equivalent  to  the  continuous  spectra,  then 


G(f,x)  = 


S(f  „x) 
(f  ,0)   k(xs9"; 


=  e 


■rax 


(17) 


which,  as  the  attenuation  coefficient,  m  =  m(f),  is  sought,  can  be 
written : 


In 


m(f4)  = 


S(f.,0)   k(x,0,f) 
S(f  .,x) 


x 


(18) 


A  valid  comparison  of  energy  densities  for  the  determination  of 
m(f .)  can  only  be  made  when  the  identical  wave  components  are  present 
at  both  measurement  points.   Therefore,  S(f  .,0)  measured  at  the  upwind 
probe  can  only  be  compared  with  S(f  .,x)  measured  at  the  downwind  probe 
after  the  time,  t,  required  for  the  energy  of  the  component  frequency 
f  .  to  propagate  over  the  distance  x,  the  effective  distance  between  probes 
This  is  given  for  gravity  waves  by 


t  -  x  /  V  =4x7Tf/g 


(19) 


where  V   is  the  group  velocity  of  the  f.  and  g  is  the  acceleration  of 
g  J 

gravity. 

Table  1  summarizes  the  data  and  computations  for  m(f  .)  for  data  set 
A.   As  the  S(f.)  were  computed  at  specific  times,  the  frequencies  whose 
propagation  times  corresponded  to  the  time  intervals  between  upwind  and 
downwind  spectra  were  computed  by  the  inverse  of  equation  (19)  and 
rounded  to  the  nearest  %  Hz  to  correspond  to  the  center  frequency  of  the 
band  represented  by  f  .,  unless  the  calculated  frequency  fell  midway  be- 
tween center  frequencies,  then  rounding  was  done  in  both  directions  and 
dual  calculations  made. 
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The  second  calculation  of  m(f .)  under  each  frequency  uses  2  S(f .,x)  in 
place  of  S(f.,x)  to  show  the  effect  of  the  possible  error  in  the  down- 
wind spectra  mentioned  in  section  A. 


TABLE  1 

SUMMARY  OF  STATISTICS  AND  ATTENUATION  COEFFICIENTS  FOR  DATA  SET  A 

t     x     f.     S(f.,0)  S(f.,x)   m(f.) 

J       J  J        J 

Upwind     Downwind                            k(x,0,f)   _        . 
time time sec    ft sec   mm  /Hz mm  /Hz   ft 

11:57:00   11:58:52   112   138   2.0   1.31      .95    .648     .00473 

-.000348 
11:57:00   11:59:56   176   138   3.0   2.24      .95    .136     .0199 

.0149 
3.5   2.35      .95    .0349    .0301 

.0251 
11:57:00   12:01:00   240   138   4.5    .519     .95    .00265   .0377 

.0327 


Now,  as  m  =  m(f),  assume  that  m(f)  can  be  expressed  by  the  equation 

m(f)  =   K  fp  (20) 

where  K  is  an  attenuation  constant  and  p  is  an  integer. 

Table  2  shows  the  values  of  K  for  the  values  of  m(f)  given  in 
Table  1  for  p  -  2,3,4. 
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TABLE  2 

m(f) 

ft-1 

ATTENUATION 

CONSTANTS  FOR 

DATA 

SET  A 

f 
-1 

sec 

m(f) 

-1 
nm 

p  -  2 

K 

2. 
sec  /nm 

p  -  3 

K 

3, 
sec  /nm 

p  =  4 
K 

sec  /nm 

2.0 

.00473 

28.4 

7.1 

3.55 

1.77 

-.000348 

-2.08 

-.52 

-.26 

-.13 

3.0 

.0199 

119. 

13.4 

4.41 

1.47 

.0149 

89.3 

9.9 

3.31 

1.10 

3.5 

.0301 

180. 

14.7 

4.18 

1.20 

.0251 

150. 

12.2 

3.49 

1.00 

4.5 

.0377 

226 

11.1 

2.46 

.545 

.0327 

196 

9.6 

2.13 

.472 

C.       COMPARISON  WITH  OCEAN  SWELL  DATA 

Some  energy  density  spectra  of  recorded  Pacific  Ocean  swell  from 
fetches  at  known  distances  from  the  swell  observation  point  were  made 
available  by  Wickham  (12).   Values  of  m  and  K  for  these  data  have  been 
computed  in  the  manner  prescribed  for  these  same  values  in  Tables  1  and 
2  above,  and  are  displayed  in  Table  3.   The  data  for  Table  3  may  include 
situations  in  which  the  winds  in  the  propagation  areas  were  across  or 
opposed  to  the  direction  of  propagation,  where  the  data  for  Tables  1 
and  2  were  for  only  the  following  wind  situation. 
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TABLE 

3 

ATTENUATION  COEFFICIENTS  AND 

CONSTANTS  FOR 

OCEAN 

WAVES 

f 

m(f) 

P  "  2 
K 

P  =  3 
K 

p  =  4 
K 

-1 

sec 

-1 
nm 

2. 
sec  /run 

3, 
sec  /nm 

4. 
sec  /nm 

.105 

.00275 

.275 

2.75 

27.5 

.095 

.00224 

.248 

2.62 

27.5 

.095 

.00261 

.289 

3.04 

32.0 

.085 

.00266 

.369 

4.34 

51.0 

.075 

.00378 

.670 

8.95 

120. 

.055 

.00115 

.380 

6.93 

125. 

D.    RESULTS  OF  DATA  SET  A  AND  TABLE  3 

1.  The  possible  error  of  a  factor  of  2  between  spectra  as  dis- 
cussed in  section  A.  does  not  have  a  significant  effect  on  the  deter- 
mination of  attenuation  for  data  set  A.  With  the  exception  of  the  values 
at  f .  =  2.0  Hz,  the  attenuation  coefficients,  m(f)  are  all  of  the  same 
order  of  magnitude;  and  with  the  exception  of  the  negative  value  for  the 
second  estimate  of  K  for  2.0  Hz,  the  attenuation  constants  for  each  power 
of  frequency  are  of  the  same  order  of  magnitude. 

2.  If  multiplying  S(f  .,  x)  by  2  had  a  small  effect,  then  variations 
within  the  confidence  limits  for  the  spectral  estimates  will  also  have 
small  effect. 

3.  It  is  quite  evident  that  if  an  attenuation  function  can  be 

described  over  a  wide  range  of  frequencies  by  a  simple  relationship  such 

3 
as  proposed  by  equation  (20),  then  it  must  be  a  function  of  f  . 
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E.  DATA  SET  B 

The  wind  in  data  set  B  died  more  slowly  than  the  dramatic  change 
which  allowed  the  extended  calculations  for  data  set  A.   For  those  short 
time  periods  where  the  wind  declined  rapidly,  S(f,)  were  computed  and 
are  displayed  in  Figures  9,  10,  and  11.   As  each  figure  represents  only 
one  attenuation  observation,  the  corresponding  values  of  m(f .)  and  K  are 
shown  on  the  figures. 

Less  confidence  is  placed  in  the  results  of  data  set  B,  as  the 
results  are  more  erratic  and  the  spectra  seem  to  be  influenced  near  the 
frequencies  of  interest  by  a  low  frequency  disturbance  that  is  not  evi- 
dent in  data  set  A.   This  low  frequency  could  have  been  due  to  a  surge 
set  up  in  the  lake  by  the  gusting  of  the  wind  during  its  decline. 

F.  ADDITIONAL  OBSERVATIONS 

The  persistence  of  the  energy  density  at  2.0  Hz  is  interesting, 

especially  in  view  of  the  negative  values  obtained  when  2  S(f  .,x)  is 

substituted  for  S(f ,,x)  for  the  calculation  of  m(f  .)  and  K.   These 

J  J 

negative  values  indicate  that  the  energy  in  that  particular  frequency 
may  be  increasing  with  distance.   A  possible  cause  for  this  is  that  2  Hz 
may  be  a  low  enough  frequency  to  suffer  negligible  attenuation  in  such 
a  short  distance. 

The  second  peak  that  occurs  in  both  data  sets  at  6  Hz  is  also  note- 
worthy.  Possible  explanations  for  this  are  that  a  frequency  of  6  Hz  is 
beginning  to  be  influenced  by  surface  tension  effects,  or  the  effect  may 
be  a  result  of  non-linear  interactions  as  have  been  reported  by  Pierson 
(7). 
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FIGURE     10.      ENERGY     DENSITY    SPECTRA -DATA    SET    B 
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FIGURE 


ENERGY     DENSITY     SPECTRA-  DATA    SET    B 
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The  decision  to  establish  the  Nyquist  frequency  at  128  Hz  was  not 
unwarranted,  as  examination  of  the  spectra  beyond  the  water  surface  wave 
frequencies  show  occasions  of  significant  energy,  particularly  around 
60  Hz. 
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VII.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  instrumentation,  recording,  filtering,  digitizing  and  analyti- 
cal techniques  employed  in  support  of  this  investigation  were  extremely 
effective  in  providing  well  defined  energy  density  spectra  from  waves  in 
a  small  scale  environment, 

The  dying  wind  situation  gives  a  definable  separation  of  wind  wave 
generation  and  swell  propagation  for  high  frequency  waves  when  the  wind 
dies  rapidly  and  smoothly, 

Strong  evidence  is  provided  in  support  of  a  contention  that  an 
empirical  attenuation  function  can  be  applied  to  a  wide  range  of  fre- 
quencies, and  it  is  proposed  that  this  function  might  be 

r.  /  r  \  — K.IX 

G(f,x)  =   e 

3 
where  K  is  approximately  3  sec  /nm  for  following  winds  and  may  increase 

3 
to  as  much  as  8  sec  /nm  for  opposing  winds.   Future  extensive  investi- 
gation may  show  the  frequency  exponent  to  not  be  an  integer,  but  for  the 
current  state  of  definition  of  the  attenuation  function  this  integral 
representation  is  satisfactory. 

It  is  recommended  that  further  studies  be  undertaken  to  provide 
better  definition  of  such  an  empirical  attenuation  function,  with  emphasis 
placed  on  intermediate  frequencies  between  the  high  frequencies  investi- 
gated in  this  paper  and  the  lower  frequencies  more  frequently  studied  in 
the  ocean. 
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